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m ; Abstract 
CN 

This work is devoted to the study of diquark correlations inside the nucleon. We analyze some 



matrix elements which encode information about the non-perturbative forces, in different color 
anti-triplet diquark channels. We suggest a lattice calculation to check the quark-diquark picture 
and clarify the role of instanton-mediated interactions. We study in detail the physical properties 



■ of the + diquark, using the Random Instanton Liquid Model. We find that instanton forces are 



sufficiently strong to form a diquark bound-state, with a mass of ~ 500 MeV, which is compatible 



Qh! with earlier estimates. We also compute its electro-magnetic form factor and find that the diquark 

<D . 

4^ ■ is a broad object, with a size comparable with that of the proton. 
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I. INTRODUCTION 



A complete description of the strong interaction between quarks is available only at 
asymptotically short distances, where the QCD becomes weakly coupled. At the pertur- 
bative level, it is well known that two quarks strongly attract each other if they are in a 
J p = + and color anti-triplet state Q, 0]. The question whether a particularly strong at- 
traction in this channel is present also at intermediate and large distances -where the theory 
is strongly coupled- has long been debated. From a phenomenological point of view, the 
notion of strongly correlated, quasi-bound diquarks has been widelv used to describe low-g 2 
processes in hadrons ( for early studies see (J, for a review see for further references 
see also lal). For example, it has been observed that the anti-triplet scalar diquarks play an 

□tin 

important role in non-leptonic weak decays of both hyperons and kaons jy, 13, |8( . 

The problem of understanding the role of diquark correlations in QCD has recently be- 
come particularly important, after several experimental groups have reported evidence for 
the first truly exotic state, the 9 + pentaquark 0|. The existence of such a state was orig- 
inally predicted by Diakonov, Petrov and Polyakov using a soliton model |h| and several 
experiments are being performed to confirm this important discovery. There are also on- 
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going theoretical inverstigations based on lattice QCD simulations (see for example 
references therein). 

The claims of the observation of the pentaquark have triggered a huge theoretical activity, 
aiming to identify possible dynamical mechanisms which could explain the existence and the 
extremely narrow width of this resonance. Jaffe and Wilczek [1Q] have suggested that several 
important properties of the hadronic spectrum, including the existence of a pentaquark anti- 
decuplet, could be explained by assuming that the non-perturbative quark-quark interaction 
is particularly attractive in the anti-triplet + channel, leading to strongly correlated scalar 
u — d diquarks, inside hadrons. 

The implications of the diquark hypothesis on the structure and the decay of the pen- 
taquark has been analyzed in a number of recent works (for an incomplete list see [3]). At 
the same time, the Jaffe- Wilczek analysis has opened a related discussion concerning the 
possible non-perturbative mechanisms in QCD, which can lead to a strong attraction in the 
scalar anti-triplet channel. Shuryak and Zahed have performed an analysis of the hadronic 
spectrum assuming that diquark correlations are induced by instantons Il3|] . In addition to 
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diquarks, Kochelev, Lee and Vento have also suggested the existence of instanton-induced 
tri-quarks Q|. It should be mentioned all models in which the quarks in the pentaquark 
are strongly correlated contrast with the original mean-field picture proposed by Diakonov, 
Petrov and Polyakov, in which two-body correlations are neglected 

Even though the experimental evidence for the + pentaquark is still a matter of debate, 
understanding if strong correlations in the + anti-triplet diquark channel are present in 
the vacuum and in hadrons is an important open problem in QCD. At present, Lattice 
Gauge Theory represents the only available framework to perform ab-initio non-perturbative 
calculations in QCD. It is therefore natural to attempt to use it to address this fundamental 
question. In this paper, we take a first step in this direction. We identify and study a set of 
lattice-calculable matrix-elements, which are very sensitive to the quark-quark interaction 
in different diquark channels. In order to show that these matrix elements can be used to 
gain insight on diquark correlations inside hadrons, we compute them using three different 
phenomenological models: a naive S77(6) Non-Relativistic Quark Model, in which dynamical 
diquark correlations are absent, the Random Instanton Liquid Model, where scalar diquarks 
are strongly favored by the 't Hooft interaction, and a Chiral Soliton Model, where the 
't Hooft interaction is treated at the mean-field level. We show that these models lead to 
very different predictions. We also identify a lattice calculation which can check if the strong 
attraction in the scalar anti-triplet diquark channel can be mediated by instantons, as in 
the model of Shuryak and Zahed. 

In the second part of this work, we study the physical properties of the diquarks induced 
by instanton-mediated forces. We provide unambiguous evidence that the 't Hooft interac- 
tion leads to a + bound-state of about 500 MeV. This result agrees with the early estimate 

n 

by Shuryak and collaborators 30J. We also compute the electro- magnetic form factor of 
such a bound-state, finding that its size is comparable with that of the proton. This result 
suggests that diquarks cannot at all be modeled as point-like objects. 

The paper is organized as follows. In section |Hl we define the relevant matrix elements 
and explain why these quantities can be used to study diquark correlations, in different 
diquark channels. The matrix elements are computed in different models in section IIHI and 
the results are analyzed in section IIVI The diquark mass and size are computed in the 
Random Instanton Liquid Model, in section All the results and the conclusions of this 
work are summarized in section IVT1 
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II. DIQUARK DENSITIES IN THE PROTON 

Much of the theoretical studies of the internal structure of hadrons focus on observables 
related to one-body local operators. These matrix elements allow to access information 
on the hadron internal wave-function, but are only indirectly sensitive to the Dirac and 
flavor structure of the non-perturbative quark-quark interaction. For example, models of 
the proton with drastically different diquark content typically give comparable charge radii. 

Much more direct insight on quark-quark correlations inside hadrons can be learned by 
focusing on two-body local density operators, which simultaneously probe the position and 
the discrete quantum numbers of two quarks. In this view, we consider a set of four-field 
operators (which we shall refer to as diquark densities), in which two quarks with quantum 
numbers of a scalar, vector, axial-vector and pseudo-scalar color anti-triplet diquark are 
destroyed and re-emitted, in the same point. 

We work with Euclidean space-time 1 , and we define: 



where C = 27274 is the charge-conjugation operator, T e {1,75,7^,7^ 75,07^} and a,b,c 
are color indexes The operator Ffi(x) absorbs two quarks in a point with quantum numbers 
of a color anti-triplet diquark. For example, F^ 3 (x) annihilates an axial- vector diquark, while 
F^ 5 (x) absorbs scalar diquarks. Hence, the matrix element: 



measures the probability amplitude to find at time t = two quarks at the point r in the 
proton, in a anti-triplet color state and with quantum numbers specified by T. 

The choice of T in determines the transformation property of the F£;(x) operator 
under exchanges of the flavor indexes. Acting on the vacuum, the operator F£ ' (r, 0) creates 
at time t — the diquark state: 



Dr(r) = F r at (r,0)F r a (r,0), 
F r a (x) = e abc u b (x)CTd c (x), 



(1) 
(2) 



Pt(t) = (P\Dt(t)\P) 



(3) 




(4) 



In our notation x 



(r,Xi). 
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A (T = 0) 
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A 



TABLE I: Symmetry properties of the operators Ffi(x), under exchange of space, color, flavor and 
Dirac indexes. 

where s, t are spinor indexes, a, b, c are color indexes, u, d are quark flavors. This state is 
symmetric under the exchange of the coordinates r and antisymmetric under the exchange 
of color. Under the exchange of the Dirac indexes, the state is either symmetric or anti- 
symmetric, depending on Y, as shown in tabled for Y = 75, the matrix Y'C is antisymmetric 
and the diquark state is thus antisymmetric under the exchange of Dirac indexes. The state 
is symmetric under the exchange of space indexes and antisymmetric under the exchange of 
color indexes. It is by construction antisymmetric with respect to the exchange of all indexes 
(space, Dirac, flavor and color) so that it is necessarily antisymmetric with respect to flavor 
indexes, meaning that it is a flavor T = state. However, when Y = 74 ( a symmetric 
matrix ) the diquark state is symmetric under the exchange of Dirac indexes. The state is 
thus necessarily symmetric with respect to flavor indexes, meaning that it is a flavor T = 1 
M T = state. 

The comparison of the different matrix elements ® can provide useful insight into the 
internal dynamics of the hadron, in particular about the strength and the spin- and flavor- 
structure of the quark-quark interaction. In general, we expect that the stronger are the 
two-body correlations inside the hadron, the larger are the discrepancies between a mean- 
field description and the exact calculation. More specifically, if the quark-quark interaction 
is particularly attractive in the scalar anti-triplet diquark channel — as it is assumed in the 
Jaffe-Wilczek picture — then we expect that the scalar diquark density p 75 (r) should be 
enhanced with respect to a model in which the interaction is not particularly attractive, in 
this channel. In fact, the u and d quarks would have a larger probability "to be found" in 
the same point and be destroyed by the local operator F* (x) . 
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The diquark densities can also be used to study relativistic effects in the hadron. The 
F±(x) and F®(x) operators mix upper and lower spinor components of the quark fields and 
vanish in the non-relativistic limit. Hence the densities pi(r) and p 74 (r) are sensitive only 
to the relativistic components of the wave-function. On the contrary the operators F® 3 (x), 
F° ' (x), and F° _{x) remain finite in the non-relativistic limit. 



III. DIQUARK DENSITIES IN DIFFERENT MODELS 

In this section, we compute the predictions for the diquark densities defined in the previ- 
ous section, using the Chiral Soliton Model, the SU(6) Non-Relativistic Quark Model, and 
the Random Instanton Liquid Model. The Non-Relativistic Quark Model is used to mimic 
a scenario in which the quarks are very weakly dynamically correlated and do not form 
bound diquarks. In such a model, all spin-flavor correlations are due to the SU(6) symme- 
try properties of the wave-function only. The Random Instanton Liquid Model calculation 
represents the opposite scenario, in which the non-perturbative interaction is spin-flavor de- 
pendent and is particularly attractive in the + anti-triplet diquark channel. The mean-field 
Chiral Soliton Model can be used to mimic a scenario in which the same interaction binds 
the nucleon, but leads to negligible two-body correlations. 

The reader who is not interested in the technical details of these calculations may skip the 
remaining part of this section and consider directly the discussion of the results presented 
in section ITVl 



Chiral Soliton Model 



The Chiral Soliton Model (for a detailed discussion we refer the reader to e.g. 21, 22]) 
can be used to account for the 't Hooft interaction in the nucleon to all orders, keeping 
only the leading terms in the 1/N C expansion. Hence, in the large N c limit, the Random 
Instanton Liquid Model and Chiral Soliton Model should give similar results. 

As all mean-field approaches, the Chiral Soliton Model is most efficient for calculations 
of matrix elements of one-body operators. Indeed, the Random Instanton Liquid Model 



and Chiral Soliton Model give similar predictions for the proton form factors 



23 



Q. On 



the other hand, if the instanton-induced two-body correlations are strong for N c = 3, then 
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we expect that the Random Instanton Liquid Model and the Chiral Soliton Model should 
give significantly different results for matrix elements of two-body operators 2 , such as the 
diquark densities (jHJ)- 

For Nf = 2, the Chiral Soliton Model is derived by bosonization of the four-field 't Hooft 
contact interaction, introducing an auxiliary chiral field. The chirally-invariant Hamiltonian 
governing the dynamics of the quarks field ip and the chiral fields a and ir a (a = 1, 2, 3) is: 



H = J d 3 r ip ] [a ■ p + gf3{a + i^ir ■ r)) ip 



, ld''v l ( Vai-..4!V^) 2 - A(^-i ^-a;:r . (5) 

where r a are isospin Pauli matrices, g and A are coupling constants and o"o is the vacuum 
expectation value of the scalar field a. In the large N c limit, one can apply the zero- loop 
approximation and treat the (a, it) fields as classical. 

The Hamiltonian (JHJ) is quadratic in the quark fields and gives rise to the single-particle 
Dirac equation 



[-i a ■ V + + z 7 5X ■ r)} |A) = e A |A), 



(6) 



where, for convenience, we have introduced dimensionless distances and fields defined as: 



x = g or r, 



a 



Xa 



CO 



(7) 



The total energy, when iV c quarks are accommodated in the J + I = valence orbital of 
energy e A , is 



E = N c ga e x + ^ J x 



(V0) 2 + (v x ) 2 + 4(0 2 + x 2 -i 

9 



(8) 



By requiring that the total energy must be stationary with respect to infinitesimal variations 
of the chiral fields one obtains the equations: 



g 2 r 
A 



+ X 2 -l) + iV c (A|x)/?(x|A) = 0, 



— V 2 Xa + 4- X a (0 2 + x 2 -l) + iV c (A|x) i(5 lh (x|A) = 0, 



9' 



(9) 
(10) 



2 For example, the Chiral Soliton Model may not be very accurate in predicting the amplitudes for weak 
decays of hadrons, which are driven by a four-field effective Hamiltonian. 
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Eq.s ©, Q an d ffTTTf) represent a self-consistent set of equations for the quark orbitals |A) 
and the classical chiral fields. The chiral soliton state is a solution of such equations, in 
which the chiral fields are assumed to be time-independent and have a hedgehog shape: 



Xa{x) = XaXW 



(11) 



In the chiral soliton model of Diakonov and Petrov j21|, the fields a and 7r a are restricted 
to the chiral circle, a 2 + tt 2 = a 2 . With the ansatz (fTTf ■ the Dirac equation (JHJ) generates a 
bound-state composed of valence and Dirac-sea quarks: 



(12) 



n at n 4 io>. 

ptdvalence \£sea 

The first product of creation operators excites three valence quarks of different colors in a 
s-state wave-function in the form: 



(xu|A) 



1 



V2 



% |x| 
G(x) 



I ) 



(X d\\) 



1 

72 



( 



|x| 



I T ) 



^a-x|T ) 



(13) 



The radial functions F(x) and G(x) are normalized according to 



(fx f F 2 (x) + G 2 (x) 



(14) 



and can be computed numerically. 

The second product of creation operators in (jllj) generates the contribution of quarks in 
the Dirac sea. In this work, we shall neglect such a contribution and retain only the valence 
part of the soliton wave-function. 

After having solved for the valence component of the soliton wave-function |$), it is 
straightforward to compute the different diquark densities: 



F r at (r, 0) F r a (r, 0)| $) , T = 1, * 7s , 74 , 7s , 7 475- 



We begin by contracting the color index a. The diquark density becomes: 



Pr (r) = (TC) SiS2 (CT) S3S4 ($ 



cf 3 1?/ 02 t?/ 2 <f 3 - (f 3 f w C2 tit 03 cf 2 



(15) 



(16) 



where si, s^, S4, S4 = 1,..,4 are Dirac indexes and 02,03 = 1,...,3 are color indexes. By 
applying Wick's theorem we get: 



^C 3 t„ C 2 ] C-2 JC 3 _ rcaf .C2j C 3 ,C2 
Ll S 1 U S2 a S3 U S4 U S1 U S2 a s 3 U S4 



= N% (A ]xsid) (xs^d |A) (// |xs 2 «) (XS3W |//) 

+ iVp (A |xsi<i) (xs 3 m |A) (/i |xs 2 w) (xs 4 <i 

— N c (A |xsirf) (XS3M |A) {fj, |xs2«) (XS4C? 

— N c (A |xsi<i) (xs 4 d |A) (fi |xs 2 m) (xs 3 m . 



(17) 



Substituting the explicit expression for the valence orbits (llHjl . and performing the summa- 
tion over the spinor indexes, it is straightforward to obtain the result (choosing x along the 
z direction): 

(18) 
(19) 



Pi( r ) 


= 0, 






P 75 ( r ) 


= N C (N C - 




(F 2 {ga T) + G 2 {ga T) 




P7 4 ( r ) 


= 2N C (N C - 


-1) 


(F(gaor) G(ga v)\ 2 
{ (^or) 2 J 


P-ft( r ) 


= N C (N C - 


D 


(F\go Q T)-G\ga,v) 


V ^(ff(r r) 2 


P7475 ( r ) 


= N C (N C - 


D 


f^g a r) -GP(ci (T r) 


V v^G^r) 2 



(20) 

) 2 , (21) 

) 2 , (22) 

where we have used the relationship between x and r given by (|7j). 

These results correspond to the different diquark densities in the soliton state, which is 
neither an eigenstate of the angular momentum nor an eigenstate of the isospin operator. In- 
stead, it is characterized by its (vanishing) eigenvalue of the grand-spin G = J + 1 operator. 
In order to make contact with the proton, one has to perform a projection onto a state with 
(T = 1/2, Mt = 1/2). The operators F^ (x) are all color zero and flavor zero operators, ex- 
cept for the operators F^ (x) and F^ (x). Thus the operators D-p (r, 0) = Fp ' (r, 0) Ff (r, 0) 
with T = (1, 75, 7475) are flavor zero operators, so that their expectation value in the nucleon 
state is equal to their expectation value in the soliton state. 

This, however, is not true for the operators F® * (x) and (x) which are flavor T = 
1, Mt = operators. Indeed, acting at time t = on a flavor zero state |0) (which is not 
the soliton state), the operator F^ produces the state 

F$ (r, 0) |0) = e abc (74 C) st |r s db;rtu c) , 

which is symmetric under the exchange of the flavor indexes. Similarly, the hermitian 
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(23) 



conjugate operator F" 4 produces the state: 

F; (r, 0) |0> = e abc (C l4 ) ts \rtdb;r sue) , (24) 

which is a flavor T — 1, My = state. Therefore the operator D 74 (r) = F^(r, 0) F^(r, 0) 
behaves, under flavor rotations as the following mixture of T = and T = 2 operators: 

(r) = -L £)T=o,m t =o (r) _ ^- D TfM T =» (r) , (25) 

where ^ and — ^/f are Clebsch-Gordon coefficients. The second term, which is a flavor T = 2 
operator, has a vanishing expectation value in the nucleon state which has flavor T = |. 
Therefore only the first term of (|23|) contributes to the diquark density in the proton: 

p 74 (r) = (P\D, 4 (v)\P) = ±={P\D T =^=\v)\P). (26) 



v/3 

,M T =0 



We need to determine the operator D^ =0 Mt=0 (r). To do this we use a simplified and more 
transparent notation: 



(x) = F${x) = F a J{x) F; (x) = = f£(x). (27) 

Then the inverse Clebsch expansion reads: 

d t=o,m t =o (r) = -i=Flt (r , 0) Fl(r, 0) - -^J(r, 0)ift(r, 0) + ^Sfr, 0) *£,(r, 0). 

(28) 

On the right hand side, the operators F^ and are the T = 1, = ±1 operators: 

= e ahc uf (x) ( 74 C) st uV (x) F a J{x) = e abc df (x) ( 74 C) st ^ (*) . (29) 

Since the ]jT=o,m t =q ^ j g a fl avor zero operator, its expectation value in the proton state 
\P) is equal to its expectation value in the soliton state |$). Therefore, the diquark density 
p 74 (r) of the proton is: 

p 74 (r) = -±=(<S>\ - ^i#(r, 0) Fl(r, 0) - ^i#(r, 0) i$(r, 0) + -L ^(r, 0) ift(r, 0)|$) 
= -i p™(r) - ~p^(r) + i p£(r). (30) 



The densities p 7 "(r) and p 74 (r) can be computed in complete analogy with the p^(r) (which 
leads to Eq. ([20))) . For a generic direction f the result is 

P- W = (r) = -2iV c (AT. - 1) F(9a ^:° T) (f. + f*) • (31) 
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Since we have chosen r along the z direction, we get a vanishing contribution. So, in the 
Chiral Soliton Model, the density p 74 (r) in the proton state (as opposed to the density in 
the soliton state) reads: 

(32) 



Notice that projecting onto the proton leads to a reduction of this density by a factor 3. 
We can repeat the same calculation to obtain the Chiral Soliton Model prediction for 



the p 73 (r) density in the proton. Also in this case, the contribution from the p"" and p\ 



,dd 

'73 



(33) 



densities vanish, if r is chosen along the z direction. The final result is 

' 3 '{ \/2( 9t r„r)2 J 

Also in this channel, the projection onto the proton has generated an extra 1/3 factor. 

In conclusion, the Chiral Soliton Model predictions for the diquark densities in the proton 
are: 



Mr 

P 7 3( r 

P7475 ( r 



= 0, 

= N C (N C -1) 



'F 2 (ga r) + G 2 (ga r) 



27V c (iV c -l) (F(ga Q T)G(ga T) 



3 V (g(r r) 2 J 

N C (N C -1) fF 2 (ga Q r)-G 2 (ga Q r) 



3 V ^(For) 2 

= N C (N C -1)( F ^^- G ^^ 
1 V2(ga r) 2 



(34) 
(35) 

(36) 

(37) 

(38) 



B. Non-Relativistic SU(6) Quark Model 



In the conventional SU(6) Non-Relativistic Quark Model, the proton wave-function is 
defined as: 



where 



I A) = ^ {xms<Pms + Xma<Pma) ® I color) ® (spatial), 



Xma = -^Ktl - IT) T), <Pma = \(ud-du)u), 



(39) 



(40) 
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Xms = ^ KTI + IT) T-2TTI ), <p MS = ±=.\( u d + du)u-2uud). (41) 

The spacial (color) wave-function is totally symmetric (anti-symmetric). 

The diquark density operators are constructed by expanding the field operators in Ff:{x) 
on a basis of the constituent quark creation/annihilation operators and taking the non- 
relativistic limit for the spinors. 

F r a (r,0) = £ u b sl (v)d c s (v)M^ s , (42) 



where 



M$ B = v],CTv s , (43) 
v s = (x-,0,0), (44) 
Xt = (1,0), xi = (0,1). (45) 

It is immediate to verify that M], 4 S = M^ s = and therefore the densities pi(r) and p 74 (r) 
vanish in the non-relativistic limit. On the other hand, one has 

Mf s = -i Xs «T2Xs (46) 
Mf s = Xs«JiXs (47) 
M]^ 5 = tXs'a 2X s. (48) 

From these relationships it follows that the operators (x) and -F° 475 (x) absorb two quarks 
when they are coupled to zero total angular momentum, while the operator F"(x) absorbs 
them in the J = 1, J z = channel. 

The matrix elements © can be computed from the expressions (}3T?j) using the anti- 
commutation relations for the quark creation/annihilation operators. Alternatively 3 , one 
can obtain them simply by taking the non-relativistic limit of the Chiral-Soliton Model 
results ( i.e. by setting G(aogr) = in Eq.s (H2H3BD ) Q- 

The result is: 

Mr) = 0, (49) 



3 We thank D.Diakonov for making this observation and pointing out an algebraic mistake in the first 
version of this manuscript. 
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p 75 (r) = N C (N C -1) J d 3 r'r(r',r,r)il;(r',r,r), (50) 

p 74 (r) = 0, (51) 

p 7475 (r) = N C (N C -1) Jd 3 r'r(v',v,r)i/;(r',v,r), (52) 

p 73 (r) = Nc{Nc o ~ 1} / d 3 rV*(r',r,r)^(r',r,r). (53) 



3 

Notice that, as long as one is interested in ratios of densities, it is not necessary to 
specify the spatial wave-function ip(ri, r 2 , r 3 ). In the simple non-relativistic model, ratios of 
diquark densities are completely determined by the color, spin and flavor structure of the 
SU(Q) nucleon wave-function. This is in general not the case in relativistic models, or in 
SU(6) breaking Non-Relativistic Quark Models, where spin and spatial degrees of freedom 
do not factorize. 



C. Random Instanton Liquid Model 

Instantons are topological gauge configurations which extremize the Euclidean Yang-Mills 
action and therefore appear as saddle points in the semi-classical approximation to QCD. 
They generate the so-called 't Hooft effective quark-quark interaction, that solves -at least 
on a qualitative level- the U(l) problem ^| and spontaneously breaks chiral symmetry [16]. 
On the other hand, present instanton models do not lead to an area-law for the Wilson loop. 

The Instanton Liquid Model assumes that the QCD vacuum is saturated by an ensemble 
of instantons and anti-instantons. The only phenomenological parameters in the model are 
the average instanton size and density: 

p = (p) ~ 1/3 fm, n = (n) ~ 1 fm~ 4 . (54) 

These value, are fixed to reproduce the global vacuuur properties ( q uarfi aud gluon coudeu- 
sates) |19|. 

The Instanton Liquid Model can be used to account numerically for the 't Hooft inter- 
action to all orders. Such calculations are performed by exploiting the analogy between the 
Euclidean generating functional and the partition function of a statistical ensemble j^j], in 
close analogy with what is usually done in lattice simulations. After the integral over the 
fermionic degrees of freedom is carried out explicitly, one computes expectation values of the 
resulting Wick contractions by performing a Monte Carlo average over the configurations 
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of an ensemble of instantons and anti-instantons. In each instanton background configura- 
tion, the quark propagators are obtained by inverting numerically the Dirac operator. In 
the Random Instanton Liquid, the density and the size of the pseudo-particles are kept 
fixed and coincide with the average values (J54j) . On the other hand, the position and the 
color orientation of each instanton and anti-instanton are generated according to a random 
distribution. 

In order to compute the diquark densities in the Random Instanton Liquid Model we 
start by considering the following Euclidean three-point correlation function: 

G r (xi,Xf,y) = (0\J a {x f )D r (y)J a (x i )\0), (55) 

where Dr(y) is the diquark density operator defined in (JTJ), and 

J a (x) = e abc u aT (x)C^d\x)u ca (x), a = l,..,4 (56) 

is an interpolating operator which excites states with the quantum numbers of the proton. 
The correlator (}55|) represents the probability amplitude to create a state with the quantum 
numbers of a proton at point Xi, to absorb and re-emit two quarks in a given diquark 
configuration at a the point y, and to finally re-absorb the three-quark state at the point Xf. 
By inserting two complete sets of eigenstates of the QCD Hamiltonian in (jo3|) we obtain: 

GT (*- »> = 5/ j£h I ■ Tr 1(01 J(X!) 1 Nip ' s ' } > 

■ (JV(p', s') I D r (y) I Af(p, s) ).< JV(P, <0I J fa) I )] + (57) 

where |iV(p, s)) denotes a proton state of momentum p and spin s and the ellipses represent 
all terms depending on its excitations (including the continuum contribution). In the limit 
of large Euclidean separations ( \xi — y\, \xi — Xf \ and \xf — y\ — ► oo), the contribution 
from the excited states to the correlation function is exponentially suppressed and only the 
proton state propagates between the operators. 

The overlap of the current operator with the nucleon can be written as: 

(0\J(x f )\N(p',s'))= A, s (p)e^, (58) 

where v s (p) denotes a Dirac spinor. Following jjjj, the matrix elements of the diquark density 
operator can be parametrized as: 

(N(p\s)\D T (y)\N(p,s)) = h r (q 2 ) v s ,(p') v s (p). (59) 
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Substituting (jB^j) and (J55|) into (joTj) one obtains: 

• /i r (g 2 ) Tr [(/ + M) + M)\ + ... (60) 

Next, we use the definition of the free fermion propagator in the forward time direction: 

S(x'; x)=f f 3 P 9 e^'-)(^ + M), (*' 4 > * 4 ) 

= 7m (x - x)^ n x (|x' — x|) + M n 2 (|x' - x|) (61) 

where 

n,(x) = ^ (k (M \x\) + ^ K X (M |x|)) (62) 
zM 2 

n 2 (x) = K\{M \x\) (63) 

and introduce the Fourier transform of the function hr{q 2 ), 

Rr(z) = j^y j*" 1 hr(q 2 ), (64) 

We obtain: 

G r (xi, Xf,y) = A 2 J d 4 z Rr(y — z) Tr [ S(xf\ z) S(z; xi) ]. (65) 

The physical interpretation of this result is the following (see also Fig. ^J. In the large 
Euclidean separation limit, the correlator G r (xi,xj,y) is governed by the function Rr(x), 
which encode information about the probability amplitude to find the diquark at a given 
distance from the center of the nucleon. Notice that Xi,Xf,z and y are four- dimensional 
vectors, so Eq. ()65j) accounts for relativistic retardation effects. On the other hand, the 
convolution of Rr{x) with the trace of proton propagators takes into account the center of 
mass motion. 

In order to clarify the relationship between the correlator (I55|) and the diquark density (jSJ) 
it is instructive to consider first the static approximation for the nucleon, in which M oo. 
In this limit, the proton propagator reads: 

S(x; y) ~ 5(x - y) e ~ M ^ &±^- t (66) 
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FIG. 1: Graphical representation of the integrand in Eq. (|65|). 

and Eq. (ffiKJ) becomes: 

G v {x i ,x f ,y) = 2 A 2 e - (xi f- x > )M 5(x f - x,) / dz 4 R r {r,z 4 ), (67) 

where we have used translational invariance to set y 4 = and we have introduced r := y — Xi, 
the distance between the center of the proton and the position where the diquark is absorbed. 
The last integral term in this expression represents the time-integrated probability amplitude 
to find a diquark at a distance r from the center of the nucleon. We can therefore identify 
this quantity with the diquark density (J3J) defined in section ITT1 

p r (r) = J dz 4 Rr{r, z 4 ). (68) 

Hence, for an infinitely heavy nucleon, the expression relating the correlation function (|55|) 
to the diquark density is simply: 

G T {x h x f ,y) = 2A 2 e-^/-^ M (5(x f -x i )pr(r). (69) 

If the nucleon mass is kept finite, there are corrections to Eq. (}6T?|) arising from replacing the 
static propagator (frJrJ|) with the exact expression (pTTj) 4 . As a result, the convolution function 
Tr[S(xf] z)S(z, Xi)] in (J65j) — which determines the position of the center of mass of the 
nucleon — is delocalized on a volume which depends on the Euclidean time r. We have 

4 We thank J.W. Negele for his observations on this point. 
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verified that for a typical mass M ~ 1 GeV and a typical Euclidean separation r ~ 1 fm 
the position of the center of mass of the nucleon is smeared on a volume of radius ~ 0.3 fm, 
centered around the origin. 

A spatial dependence of the point-to-point propagator is a signal that our nucleon is 
not at rest. Indeed, in order to isolate completely the zero-momentum component of the 
nucleon wave-function in our calculation, one would need to use very large Euclidean time 
separations. However, this would be extremely computationally demanding, because the 
signal-to-noise ratio for the propagators drop exponentially with the Euclidean time sepa- 
ration. 

From these considerations it follows that one should not quantitatively compare the de- 
tailed shape of the diquark densities of the Chiral Soliton Model and Constituent Quark 
Model with that of the Random Instanton Liquid Model. However we shall see below that, 
as long as one is interested in ratios of diquark densities, it is still possible to draw important 
qualitative conclusions about the relative strentgh of the diquark correlations. 

In this work, we shall estimate ratios of different diquark densities pr'( r )/Pr(r) by com- 
puting ratios of correlation functions G r ' {xf 1 Xi,y)/G T {xf 1 Xi 1 y). Taking ratios allows to 
maximally reduce the effects of the corrections due to center of mass motion. In fact, 
such corrections affect in the same way both the numerator and the denominator and do 
not change the normalization of the correlation function. We have computed choos- 
ing Xj = Xf = 0, y = (0,0, r, 0) and x^ = —x\ := r, with r = 1 fm. Previous anal- 
ysis in the Random Instanton Liquid Model have shown that the proton pole is isolated 
from its excited states for r ^ 0.9 fm 

0, H Q- 

The numerical calculation has been 
done averaging over 800 configurations of an ensemble of 492 instantons of 0.33 fm size, in 
V = 4.5 3 x 5.4 fm 4 box with the topology of a torus. We have used rather large current 
quark masses (m u = rrid — lOOMeV) in order to reduce finite volume artifacts. 

Let us comment on the dependence of the Random Instanton Liquid Model predictions 
on the value of the bare (current) quark mass chosen in our simulations. It is well known 
that quarks propagating in the instanton vacuum acquire an effective mass associated to 
the breaking of chiral symmetry. More specifically, the position of the pole of the quark 
propagator in the instanton background is shifted by an amount M e ff(p, m q ), which depends 
on both the momentum of the quark p and of its current mass, m q . 

The dependence of the effective mass on the current quark mass m q has been studied in 
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a number of works, using different approximations (see e.g. |35j and references therein). In 
all such studies it was found that the sum of the effective mass at zero momentum and the 
current quark mass M const = m q + M e ff(0, m q ) (which can be interpreted as the constituent 
quark mass) remains of the order ~ 350 MeV for all current masses m q < 200 MeV. In 
other words, in this model the constituent quark mass is rather insensitive on the the value 
of the current quark mass. So, we do not expect that the predictions of Random Instanton 
Liquid Model presented in this work would change dramatically if a smaller current quark 
mass was used. 



D. Note on Lattice Calculation of Diquark Densities 

We conclude this section by noting that the calculation scheme based on point-to-point 
correlation functions, used to compute the diquark densities in the Random Instanton Liquid 
Model, can also be applied to compute the same quantities in Lattice QCD. Conceptually, 
one needs to replace the average over the configurations of the instanton ensemble with an 
average over a/Hattice configurations, performed in the usual way, i.e. by sampling the space 
of lattice links. Unlike instanton models, lattice QCD calculations are affected by ultraviolet 
divergences, which are regularized by the lattice spacing. Hence, in a lattice computation, 
the diquark density operators -D r (r) have to be treated as usual Wilson operators and need 
to be renormalized. 



IV. COMPARISON AND DISCUSSION 



In this session, we discuss the results of the phenomenological calculations presented in 
the previous session. We choose to consider only the ratios constructed by dividing the 
densities p 73 (r). p 74 75 (r) and p 74 (r) by the scalar density p 75 (r). There are several reasons 
for this choice. In general, ratios are less model dependent than the individual densities. For 
example, in the Non-Relativistic Quark Model, they are completely insensitive to the details 
of the spatial wave function. Moreover, in lattice and Random Instanton Liquid Model 
calculations, taking ratios allow to remove the dependence on the Proton Mass, Euclidean 
time and the coupling A (see Eq. (JHHJ)). We shall not discuss the ratio constructed with the 
pi (r) density, as it is identically zero in all the models we have considered. 
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FIG. 2: Predictions for the ratio i?i(r) = /0 73 (r)/p^ 5 (r) in the Non-Relativistic Quark Model 
(NRQM), in the Random Instanton Liquid Model (RILM) and in the Chiral Soliton Model (CSM). 
The CSM results include the projection onto the proton state. 

Let us first analyze the ratio constructed with the density of axial-vector and scalar 
diquarks: 

= (70) 

The results of our phenomenological calculations are shown in Fig. El In the SU (6) Non- 
Relativistic Quark Model, this ratio is completely determined by the SU(6) spin-flavor 
structure of the wave-function, and is identically equal to 1/3. In the Random Instanton 
Liquid Model, -Ri(r) is sizebly reduced in magnitude (by a factor ~ 5). In the Chiral 
Soliton Model, at r = the p-wave contribution from the lower components of the spinors 
vanishes and one recovers the Non-Relativistic Quark Model Results. On the other hand, 
the axial-vector diquark density drops down very rapidly at the border of the soliton, where 
the pion field is most intense. 

These results can be interpreted as follows. In the Random Instanton Liquid Model, 
the spin- and flavor- dependent 't Hooft interaction generates a strong attraction which 
enhances the probability amplitude of finding two quarks in the same point in the + anti- 
triplet configuration, relative to the amplitude of finding them in the 1 + configuration. This 
explains why the Random Instanton Liquid Model prediction for Ri(v) is much smaller 
than that of the 577(6) Non-Relativistic Quark Model and the Chiral Soliton Model. The 
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discrepancy between Random Instanton Liquid Model and mean-field Chiral Soliton Model 
calculations is a signal that in the instanton vacuum at iV c = 3 there are strong two-body 
correlations, which are not captured by the large N c approximation. 

From this comparison it follows that a lattice calculation of i?i (r) could provide informa- 
tion about the strength of scalar diquark correlations in the nucleon. If the non-perturbative 
QCD interactions generate a strong correlation in the + anti-triplet channel, as assumed in 
the Jaffe-Wilczek model, then we predict that the curve obtained from a lattice calculation 
should lie much below 1/3. If lattice simulations found that R\{r) ~ 1/3, then this would 
imply that diquarks are not particularly correlated in the + diquark channel and therefore 
the Jaffe-Wilczek picture is not correct. -Ri(r) > 1/3 would represent an indication that the 
quark-quark interaction is less attractive in the + channel, relative to the 1 + channel. This 
would certainly be a very surprising result, since the (3 C , 3/) channel is favored by both the 
perturbative and the instanton-mediated interactions. 

Unfortunately, the ratio -Ri(r) does not encode information about the microscopic dy- 
namical mechanism underlying such diquark correlations. In fact, two completely different 
quark-quark effective interactions (e.g. one with a chirality-conserving vertex and one with 
a chirality-flipping vertex) may lead to the same predictions, as long as the short-range 
attraction in the scalar channel is sufficiently strong. 

In order to gain some insight on the microscopic origin of diquarks we need to analyze a 
different ratio: 



The results of our calculations in the three phenomenological models are reported in 
Fig. El In the Non-Relativistic Quark Model, both the p 75 (r) and p 75 74 (r) densities probe 
the + scalar diquark content of the proton, so -R2(r) = 1. In the Random Instanton Liquid 
Model the magnitude of this ratio is smaller than in the Non-Relativistic Quark Model, by 
a factor 3 or so. The fact that, in the Random Instanton Liquid Model, -R2( r ) ^ 1 has 
an important dynamical explanation. It is due to the different sensitivity of the numerator 
and denominator to the so-called direct-instanton contribution. The p 75 (r) diquark density 
receives maximal contribution from the interaction of quarks with the field of the closest 
(direct) instanton, in the vacuum. This statement can be verified by computing the cor- 
relator in the single- instanton approximation, discussed in j^J. On the other hand, the 




(71) 
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FIG. 3: Predictions for the ratio i?2(r) = /0 7475 (r)//9 75 (r), in the Non-Relativistic Quark Model 
(NRQM), in the Random Instanton Liquid Model (RILM) and in the Chiral Soliton Model (CSM). 



density in the numerator, p 75 74 (r) does not receive such a direct- instanton contribution and 
instanton-induced effects come only from the interactions of quarks with many instantons. 
The magnitude of the latter contributions are parametrically suppressed by the diluteness 
of the instanton vacuum, k = np 4 <C 1. Physically, this is the same reason why vector and 
axial- vector channels have a rather large "Zweig rule", forbidding flavor mixing, while for 
scalar and pseudo-scalar channels such a mixing is very strong. 

We remark that the very strong channel dependence of hadronic correlation functions is a 
well-known dynamical implication of instanton models. It is quite hard to obtain this effect 
in alternative dynamical mechanisms (for a detailed discussion see j26|). This was initially 
pointed out by Novikov, Shifman, Vainstein and Zakharov, in the context of O per ator Prod- 
uct Expansion (OPE) and QCD sum-rules for hadronic correlation functions 27[. Diquark 
correlations induced by the direct instantons were discussed also in j^ . 

In the Chiral Soliton Model, the ratio -R2( r ) remains of order 1 for |r| < 1 fm and 
drops rapidly at the border of the soliton. The significant deviation of the Chiral Soliton 
Model result from the Random Instanton Liquid Model prediction shows that a mean-field 
approach does not capture correlations associated to the direct-instanton effects. 

From this discussion it follows that, if the scalar diquark correlations are mainly in- 
duced by instantons or in general by a Nambu-Jona-Lasinio type of interaction (i.e. chi- 
rally symmetric, with a chirality flipping vertex), then a lattice measurement should give 
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FIG. 4: Predictions for the ratio i?3(r) = p 74 (r)//3 75 (r) in the Non-Relativistic Quark Model 
(NRQM), in the Random Instanton Liquid Model (RILM) and in the Chiral Soliton Model (CSM). 
The CSM results include the projection onto the proton state. The Non-Relativistic Quark Model 
curve is identically zero. 



The results of our calculations in the three phenomenological models are reported in Fig. 0] 
This quantity is identically zero in the Non-Relativistic Quark Model, so it represents a 
probe of relativistic effects in the nucleon. The ratio is also rather small in the Random 
Instanton Liquid Model, denoting that quarks propagating in the nucleon receive small 
relativistic corrections 5 . On the contrary, Rs(r) in the Chiral Soliton Model increases rapidly 
and approaches 1/3, near the border of the soliton. The interpretation of this result is 
the following 6 . Near the edge of the soliton, the quark field changes very rapidly. This 
rapid change gives rise to a large derivative of the wave-function in coordinate space or, 
equivalently, to large high-p components, in momentum space. Such modes enhance the 
contribution from the lower components of the spinor wave-function, which give rise to 

5 We recall that in both the Random Instanton Liquid Model and in the Chiral Soliton Model relativistic 
effects are included. 

6 We thank W. Weise for his comments on this point. 



p 7475 (r) < p 75 (r), so i? 2 (r) < 1. 
Let us now analyze a third ratio: 




(72) 
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relativistic corrections. 



V. PHYSICAL PROPERTIES OF THE SCALAR DIQUARK IN THE INSTAN- 
TON VACUUM 

In the previous section we have shown that instantons generate a strong attraction be- 
tween a u and a d quark in the color anti-triplet + channel. In this section, we provide 
unambiguous evidence that such forces are strong enough to form a bound diquark state. 
We also estimate the diquark size by computing its electric charge radius in the Random 
Instanton Liquid Model. 

A. Mass of the scalar diquark 

The question if instanton models lead to a bound + color anti-triplet diquark has been 
first posed by Diakonov and Petrov j3] and investigated in a number of works. In an 
exploratory study j^J Shuryak, Schafer and Verbaarshot computed some diquark two-point 
functions in coordinate space, using the Random Instanton Liquid Model. They found some 
evidence for a light diquark bound-state (of rougly 450 MeV mass) , by performing an analysis 
of the correlation function, based on a pole-plus-continuum parametrization of the spectral 
dens.ty. On the other hand, D.akonov, Petrov and collaborators have ana,yzed diq „ar ks by 
solving Schwinger-Dyson equations at the leading order in 1/N C [31]. They found evidence 
for correlations, but no binding 7 . 

In order to clarify this issue, we have followed an approach which is usually applied to 
extract hadron masses in lattice simulations. As usual, we have replaced the average over 
all gauge configurations with an average over the configurations of the instanton ensem- 
ble. This method presents some advantages with respect to the coordinate-representation 
calculation of Shuryak et al. On the one hand, it avoids the undesired additional model de- 
pendence associated with the parametrization of the spectral function. On the other hand, 
it makes unambiguously evident the existence of the bound state and allows to determine 
more precisely its mass. 

7 Except for N c = 2, in which case the diquark is a baryon and its mass is protected by Pauli-Gursey 
symmetry. 
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The starting point consists of computing the Euclidean two-point function: 




(73) 



where J_d(x) is the usual diquark current: 



J' 



( x ) = e abc u bT (x)C-f 5 d c (x). 



(74) 



The path-ordered exponent in (J73J) represents a Wilson line connecting the two extremes 
of the two-point function and is needed to assure gauge invariance of the correlator. In 
the instanton vacuum, the contribution from the Wilson line is very small, as heavy quarks 
couple very weakly with instantons. 

In the large Euclidean time limit, only the lightest state with the quantum numbers of 
the current Jd{%) propagates in the two-point function. If there is a diquark bound-state, 
then, in the large Euclidean time limit, the logarithm of the two-point function must scale 
linearly with r: 



where is the mass of the diquark and the constant is its coupling to the current, 
defined by (0\J D (0)\D) = A D . 

We have computed such a correlation function in the Random Instanton Liquid Model, 
averaging over 100 configurations in a 4 3 x 8 fm 4 box. In analogy with lattice simulations, 
a rather large quark mass was needed to reduce the finite- volume artifacts 8 . The result for 
our calculation of the correlation function ()73|) with m q ~ 100 MeV is shown in Fig. |SJ 
We clearly see that, in the large r limit, the logarithm of the two-point function becomes a 
linear function of r. 

In order to study the sensitivity of the diquark mass on the current quark mass, we 
have computed the correlation function for different values of the u and d quark mass, 
m u = md — m q — {180, 140, 100} MeV. The dependence of diquarks mass on the current 
quark masses is plotted in Fig. |H1 A rough estimate of the diquark mass can be obtained by 
linearly extrapolating to the physical value of the current quark mass 9 . We find a diquark 
mass of ~ 500 MeV, in good agreement with the previous estimate. 

8 We choose quark masses of the same order of magnitude of the smallest bare masses used in present lattice 
simulations. 

9 A precise extrapolation fit should include logarithmic corrections associated with chiral physics. We do 
not need to discuss these effects here, as we are interested in the order of magnitude of the diquark mass. 



lnG2(r) ™lnA^-r M< 



D, 



(75) 
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FIG. 5: Logarithm of the diquark two-point function, lnG2(r) computed in the Random Instanton 
Liquid Model. The linear slope is a clean signature of the existence of a bound state. 

B. Diquark size 

An important question to address is the size of the diquark. From general arguments, 
we expect that it should be comparable with the size of the proton. In fact even the most 
tightly bound QCD excitation, the pion, is known to have a rather large electric charge square 
radius, r£ ~ 0.6 fm. On the other hand, an early naive analysis in the Instanton Liquid 
Model 3^ gave a surprisingly small estimate for the diquark charge radius, te ~ 0.3 fm. 
It is hard to estimate the accuracy of such an estimate which was based on a number of 
assumptions 10 . Hence, it is worth performing a direct calculation of the diquark electro- 
magnetic form factor. 

To this end, we have computed the diquark electro-magnetic three-point function, 

G3 D (r,q) = J d 3 yd 3 x e l (0| T[ J D (0,r) J 4 e/m (y,0) f D (x,-r) Pel ]|0) (76) 

where Jd(x) is the usual diquark interpolating operator and Jf/ m {x) =2/3 u{x) 7 M u{x) — 
1/3 d(x) 7 M d(x) is the electro-magnetic current operator. In the large r limit, the correlation 
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For example, it was assumed that the diquark form factor follows a monopole fit 
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FIG. 6: The mass of the u-d diquark, computed in the Random Instanton Liquid Model, as a 
function of the u and d bare quark mass. 



functions (fTrjj) relates directly to the diquark form factor 11 , F D {Q 2 ): 

G3 D (r,q) A 2 D e -(M+V^i> Fp(Q 2 ), 



(77) 



where the values of the diquark mass and the constant A^> can be extracted from the two- 
point correlation function. 

The result of the Random Instanton Liquid Model calculation of the form factor (nor- 
malized to the total u — d diquark charge, eo = 1/3 ) is presented in Fig. where it is 
compared to the phenomenological dipole fit which reproduces the low-energy data on the 
proton electric form factor. These results imply that the diquark is an extended object, 
whose size is comparable with that of the nucleon. Its electric charge radius is found to 
be r|n (diquark) ~ (0.70 fm) 2 , which should be compared to the value of the electric charge 
radius of the proton, r|(proton) ~ (0.76 fm) 2 computed in the Random Instanton Liquid 
Model in S^j]. As expected, the electric charge radius of the diquark is larger than that of 
the pion and smaller than that of the proton. 



11 For a discussion of these relationships and of the details of the numerical integration method used for 
performing the Fourier Transform in H76[) see [2^, [3^] . 
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VI. CONCLUSIONS 



In this work we have addressed the question of how is it possible to study two-body 
diquark correlations in hadrons, using lattice QCD. We have identified some suitable lattice- 
calculable correlation functions, which allow to probe directly the diquark content of the 
nucleon. We have argued that these Green's functions provide important information about 
non-perturbative correlations, in different color anti-triplet diquark channels. 

In particular, the ratio Ri(r), defined in Eq. (|7U|h measures the strength of the correla- 
tions in the scalar diquark channel, relative to that in the axial-vector channel. It can be 
used to check the main dynamical assumption of the Jaffe-Wilczek model. In fact, if the 
quark-diquark picture is correct, then we predict that a lattice measurement must lead to 
i?i(r) < 1/3. 

The ratio -R2( r ), defined in Eq. (J7TJ), can be used to gain some insight on the dynamical 
origin of the non-perturbative interaction in the scalar diquark channel. In particular, it can 
be used to check the hypothesis according to which these forces are mediated by instantons 
(or more generally by a NJL-like, chirality-flipping interaction). We have argued that, in 
this case, we expect that lattice measurements should give i?2( r ) "C 1. 

We have computed these ratios using three phenomenological models. We have found 
that they lead to radically different predictions for the matrix elements we have selected. 
Hence, a lattice measurement could point out which picture is most realistic. It is worth 
stressing that from the fact that the Chiral Soliton Model predictions differ substantially 
from the Random Instanton Liquid Model results one should not conclude that the Chiral 
Soliton Model is not correctly reproducing the physics of instanton-induced interactions. In 
fact, such mean-field model can account for one-body local operators. On the other hand, 
our results show that it is much less reliable for computing matrix elements of two-body 
operators. The discrepancy between Random Instanton Liquid Model and Chiral Soliton 
Model should be taken as an indication that rather strong direct two-body correlations are 
generated in the instanton vacuum. 

We would like to stress that the qualitative differences between the predictions of the 
three models considered in the present analysis are dramatic and therefore quite robust. 
However, one should be careful in comparing the details of the curves associated to the 
diquark density ratios, at a quantitative level. In fact, on the one hand, the Random 
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Instanton Liquid Model predictions are affected by the smearing effects associated to the 
center of mass motion (see the discussion in section IlII Cj) . On the other hand, the Chiral 
Soliton Model and the Random Instanton Liquid Model calculations are performed using 
different values 12 of the current quark mass m q . This fact may reflect itself in a slightly 
different value of the effective mass of the quarks propagating in the chirally broken vacuum 
(see the discussion at the end of section IIII Cj) . 

In the second part of the work, we have studied in detail the physical properties of 
instanton-induced diquarks. We have provided unambiguous evidence that instantons gen- 
erate a scalar diquark bound-state of mass ~ 500 MeV, in good agreement with earlier 
estimates based on point-to-point correlation functions in coordinate space. 

We have also studied its electric charge distribution. Our results show that the scalar 
diquark is an extended object, whose size is of the order of the fm, hence comparable with 
that of the proton. Thus, phenomenological quark-diquark models cannot treat the diquark 
as a point-like object. 

As a concluding remark, we stress that the analysis performed in this work focused on 
the diquark content of the proton. However, the same study could be repeated, in order 
to investigate the diquark content of other lowest-lying baryons. For example, it would be 
particularly interesting to compare the diquark densities in the nucleon and in the delta. In 
fact, in the delta, the leading direct-instanton effects are suppressed, due to the Dirac- and 
flavor- structure of the 't Hooft interaction. Hence, if diquark correlations are instanton- 
mediated, then we expect that the strong enhancement of the scalar density observed in the 
proton should be much less pronounced in the delta. 
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FIG. 7: The diquark form factor (normalized to the total charge) evaluated in the Random Instan- 
ton Liquid Model (points) and compared with the dipole parametrization of the low-energy data 
on proton electric form factor: F^ip = 1/(1 + Q 2 /0.71) 2 (dashed line). The numerical simulation 
was performed in a (4 3 x 6) fm 4 box and the exchanged momentum is quantized accordingly. 
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